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Improved Multi-Channel Grating Design Technique 

Field of the invention 

The present invention relates broadly to a method of improving a grating design function 
describing a refractive index variation defining a multi-channel grating structure in a waveguide 
matCTial, to a metihod of fabricating a multi-channel grating structure, and to a multi-channel 
gratiAg structure. 

Badcground of the invention 

Mxilti-chaimel grating structures are typically written into photosensitive waveguides. 
The grating structure comprises a refractive index profile induced in the photosensitive 
waveguide, which in turn determines the optical characteristics such as the reflection, 
transmission, and group delay characteristics of the resulting grating structure. 

The amplitude and phase of the refractive index profile can be described by a grating 
. . design function, which in turn is used in the control of an apparatus for writing the grating 
structure into the photosensitive waveguide. Typically, the writing apparatus comprises an 
optical interferometer. 

It is known that a maximum refi:active mdex contrast required to write multi-channel 
gratings increases witti mcreasing number N of channels. Since photos^sitive materials, such 
as a photosensitive fibre used to &bricate a Bragg grating, have material limits of the maximum 
achievable photo-induced refractive index change An, this represents a limitation on the 
maximum number of channels N that can be recorded in a given material. Thus, it is highly 
desirable to reduce a required Aw^v as much as possible. 

At least preferred embodiments of the present invention seek to provide a method of 
improving a grating design fimction in terms of a reduced maximum refractive index variation 
reqmred in tiie waveguide material along the gratmg stracture, without compromising the 
quality of spectral characteristics. 

Summary of the invention 

In accordance wilh a first aspect of the present invention there is provided a method of 
improvmg a grating design fimction describing a refractive index variation defining a multi- 
channel grating stmcture in a waveguide material, the improvement being a reduced maximum 
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refiactive index variation in the waveguide ' material along the grating structure while 
maintaining a desired functional spectral domain in a spectral response function associated with 
tile design function, the method comprising the steps of: 

- modifying a first design function to generate a second design function having a reduced 
maximum ampUtude compared witii flie first design fimction, 

- detennining a second response function associated vdthihe second design fimction, 

- modifying the second response fimction to create a tiiird response fimction having a 
desired functional spectral domain, and 

" detennining a tiiird design fimction associated wifli flie fliird response fimction, 

and iterating flie method steps until tiie desired improvement is achieved. Wherein tiie 

tiiird design fimction of tiie previous iteration takes tiie place of tiie first design fimction of tiie 

next 

Jn one embodiment, tiie step of modifying tiie second i^onse fimction comprises 
rqilacing tiie corresponding spectral domain of tiie second response fimction by tiie desired 
functional spectral domain. 

The desired fimctional spectral domain may conqjrise a cort^spondin^ spectral domain 
of a first response fimction associated witii tiie first design fimction. 

The metiiod may comprise flie pre-iteration step of determining tiie first design fimction 
6om a or the associated first re^onse function. 

Preferably, flie determining of tiie response fimctions fiom tiie design fimctions 
comprises solving a direct scattering problem, and tiie detemiining of tiie design fimctions from 
tiie response fimctions comprises solving an inverse scattering problem. 

In one embodiment,, tiie tiiird response fimctiotf fi^tiier has a desired response 
characteristic in at least one fiirther spectral domain otiier tiian tiie fimctional spectial domam. 

Preferably, tiie step of modifying tixe first design fimction comprises applying a 
nonnalisation process or an averaging process to flie first design fimction. 

In one embodiment, tiie normalising process comprises replacing an ampUtude fimction 
H^) of tiie design fimction by aproduct of tiie square root of a constants and a corresponding 
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single-channel seeding amplitude function a:,(z), while maintaining a phase function of the 
design function. The constant A may be defined by the nomialisation condition 

A- ^K^dzl jK^dz , where / is a length of the multi-channel grating structure. 

0 0 

In one embodiment, the averaging process comprises averaging over a sampling or 
quasi-sampling period of the design function. 

In accordance with a second aspect of the present ifivention, there is provided a method 
of fabricating a multi-chaimel grating stmcture, the method conq)rising the step of improving a 
grating design function describing a refractive index variation defining the multi-channel 
grating structure in a waveguide material as defined in the first aspect 

In accordance with a third aspect of the present invention, there is provided a multi-* 
channel grating stmcture &bricated utilismg a method of fabrication as defined m the second 
aspect 

Brief description of the drawings 

Preferred forms of the present invention will now be described with reference to the 
accompanying drawings. 

Figures 1 (aHd) show calculated the first spectral characteristics and the first design 
amplitude of a multi-channel grating stmcture (the first design phase is not shown) 

Figures 2 (a)-(d) show calculated the second spectral characteristics and the second 
design results (normalised/averaged amplitude profile) 

Figures 3 (a)-(d) show calculated the third spectral characteristics and the third design 
results which correspond to the first completed iteration embodying the present invention. 

Figures 4 (a)-(d) show calculated spectral characteristics and design results after. 8 
iterations of an improvement process embodying the present invention. 

Figures 5 (a)-(d) show calculated spectral characteristics and design results after 150 
iterations of an improvement process embodying the present invention. ' 

Figure 6 is a flowchart illustrating a multi-chaimel grating design improvement process 
embodying the present invention. 
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Figure 7 shows an experimental set-up for writing a multi-channel grating structure of a 
mxilti-chaimel grating design embodying the present invention. 

Detailed description of the embodiments 

The preferred embodiment described provides a method of improving a grating design 
function in terms of a reduced maximum refractive index variation required in a waveguide 
material along the grating stmcture. 

Figures l(aHd) show calculated spectral-characteristics and design of a 17-channel 
Bragg grating prior to iteration procedure. It will be appreciated by a person skilled in the art 
that the spectral-characteristics are of a superb quality with virtually no deviations from a 
desired spectral-characteristics, in the example shown in Figure 1 a square-shaped transmission 
characteristic (compare Figure lb). 

It has been found by the inventors that the design of the 17-channel Bragg grating shown 
in Figure 1 can be inqproved in terms of minimisation of a maximum refractive index change 
required in fhc photosensitive material, e.g. a photosaositive optical fibre, into which the 17- 
. channel Bragg grating is to be written. In other words, the fibre Bragg grating (FBG) ampUtude 
versus distance z along the FBG (Figure la), can be improved in terms of ntiinimisation of the 
maximum amplitude without compromising the spectral response. In the example embodiment, 
this is achieved by subjecting the design to an iterative optimisation procedure which will now 
be described. 

In the example embodiment, the maximum in the FBG antrplitude versus distance z is 
initially reduced in the following way: Starting fiom the FBG amplitude versus distance z 
shown in Figure 1(a), the amplitude function k{z) is replaced by ^fAK^{z)y where ic, (z) is the 
corresponding single channel (seeding) grating amplitude function and the constant A is defined 

by the normalisation condition A= ^K^dz/ ^K^dz . Alternatively k{z) is replaced by its 

0 0 

averaged (smoothen) version. The unchanged multi-channel grating phase and the modified (the 
second) grating amplitude are then used as the input data for solving a direct scattering problem 
utilising a suitable direct scattering solver algorithm. Thus, a spectral response function 
corresponding to the modified grating design function is determined. 
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The second grating design amplitude and the coiresp<»iding spectral response data are 
shown in Fig. 2(a-d). Clearly spectral characteristics in the central part of wavelength range, 
which is the functional spectral domain in the example embodiment, are less than perfect As a 
next step these characteristics are modified. 

Within the central channel spectral domain it is replaced by the initial ^ire-iteration) 
data, compare Figure 1(b) and 3(b)- Outside the central channel spectral domains the 
deteraiined spectral response data is left intact compare Figure 2(b) and 3(b), i.e. that spectral 
domain is allowed to evolve during the iteration process. 

As a final step of each iteration, an inverse scattering problem is solved for ttie modified 
spectral response data to determine the associated, better optimised, multi-channelled grating 
design function, see Fig. 3(a-dr) 

The above method steps are then iterated in the example embodiment, until a desired 
improvemoit is achieved. Figures 4 and 5 show the level of convergence after 8 and 150 
iterations for the 17-channel grating design. For each of the Figures 4 and 5, (a to d) show 
calculated spectral characteristics and design of the 17-channel FBG. As can be seen from a 
comparison of Figure 1(a) and Figure 5(e). the optimisation process ranbodying the present 
invention has improved by about 13.5% to a lower vahie, while, due to flie nature of the 
optimisation process embodying the present invention, the superb spectral characteristics of the 
design (compare Figures 1(b), (d), and Figures 5(b), (d)). 

Figure 6 shows a flowchart illustrating the iterative process of the example embodiment, 
as follows: 

Step 500, applying a normalisation or averaging process to a givai FBG amplitude 
function; 

Step 502, solve a direct scattering problem based on the modified grating design 
function; 

Step 504, modify the resulting spectral response; 

Step 506, solve an inverse scattering problem based on the modified spectral response; 

Step 508, determine whetiier a desired quality has been achieved, and 

repeat the steps (retum loop 510) if not, or stop the iteration as indicated by arrow 512. 
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In the foUowing, a description wiU be given of how to obtain a suitable initial (pre- 
iteration) grating desiga (compare Figure 1), in example embodiments. 

In the example embodiment, a multi-channel grating function is constructed by solving 
the standard inveise scattering problem for muW-channel wavelengfli-shifted spectral 
characteristics. Spectral response functions (X) for partial single gratings are being de- 
phased with respect to each other. In other woids. the inveise scattering problan may be 
presented as: 

with nonzero relative phases 5, in the example embodiment. It is to be noted that 
spectral responses are being de-phased. not partial seeding gratmgs themselves. After solving 
the inverse scattering problem for ">"" (Pl). the multi-chamiel gratmg design fonction is 
obtained and may be presented in a form: 

M^ere we e^Ucitly retain, for illustrative purposes, the single channel grating design 
function K.e^'^^K The remainmg fectois in the expression (2) represent a "sampling" 
fimction, which is periodic. It is thus not a sampling, i.e. periodic, fixnction anymore, but 
deviates from a sampling function. 

It will be appreciated by a person skiUed in the art that in order to initially (prior to our 
iterative scheme) pre-optimise the gratmg design for minimum refractive mdex change as 

a function of the number of channels suitable numerical or analytical methods can be 
applied. 

For relatively small number of channels N one can numerically scan flirough possible 
combinations of relative phases ^ (solving an inverse scattering problem for each particular 
combmation of dqjhasmg angles Qi) and selecting the combination which is optimal according 
to some specific selection criterion (e.g. selecting the combination which minimises maximimi 
-required refractive index change). 

For N»\ location of tiie optimal set &i is difiBcult. Even rough direct scanning 
through all possible sets of angles (foUowed by efficient nmnerical mmimmn search routines) 
quickly becomes numerically mefGcient. 
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In the example embodiment, the approximate equivalence between partial spectra 
dephasing angles i9/ and partial grating relative phases ^/ is utilised. Indeed, for weak gratings 
the first order Bom approximation holds: 

-|^(z/2)=£K)9)expHAm (3) 

where q{z) is a grating design function and r{fi) is a complex reflection coefficient. 
The Fourier transform (3) is a linear operation with a major property 

F(aj /^^ + 02 r^^) ^ ai F 1^^^ ^ a2 F i^K Thus, in this approximation, dephasing of partial 
gratings is equivalent to dephasing of partial spectral channel responses. Formally the last 
statement does not hold beyond weak grating limit. However, in practice, it is still 
approximately correct and the optimal set of angles ^ (for dq>hasing of partial gratings) may be 
used as a very good approximation for the optimal set of partial spectral channel angles The 
only important change related to the ^'sampling*' function between dqphasing of partial channels 
and dephasing of partial gratings approaches is that it becomes sligjhitly. aperiodic, as mentioned 
above. 

Therefore, the phase shift values i9/ may be taken firom a partial gratings de-phasing 
grating design method In the example embodiment, a sampling function which periodically 
modulates the amplitude of a given single-channel grating (seeding gratiag) is utilised In 
addition to the periodic modulation of the amplitude of the seeding grating, different relevant 
phases ^ for each of the wavelength-seeding gratings are introduced. Accordingly, the 
resulting design function in the example embodiment may be expressed as: 

(4) 

.where the additional phase of the grating ^/^y/{z) and the sampling amplitude Q-Q(z) 
are given by: 

N/2 . . 

g^(z) = 42cos(a; -a^jcos(n;AAz/2 + P/)cos(/2^AAz/2 + p^), and 
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tff(z) = tan"' 



^^^^ sina,cos(/i;Afe/2 + y3/) 



^^^j cosa;Cos(«,Afe/2 + >S,j 



is even 



or 



(Ar-l)/2 

Q^(z) = 4 ^cosa^cos(ii;Afe/2 + p|) + 



(JV-l)/2 



4 ^cos(a, ~a^)cos(/i;Afe/2 + P;)cos(«^Afe/2 + p^) + l, 



and 



v|/(z) = tan'* 



^ J sin cos(«/ Afe / 2 + ) 



Nis odd. 



S/^r^^'^^^^i cos(«,Afe/2 + P; +1) 
where n, s2/ — iV-l and s2/?-i\^-l. 

In the above expressions for Q(z) and f//(z) we use notations 
a 5 + ) / 2, /?/ s (^^ - ) / 2 and set ^(^+,)/2 = 0 for odd number of channels. 

In order to reduce flie maxinoiiim value of Q along the grating structure, while trying to 
avoid touching the zero level along the gratmg structure, the minimization strategy used in the 
example embodiment is minimizing its maximum deviations of Q(z) along z from the 

theoretical limit level of yfW . Mathematically this may be formulated as finding 

Q (z'a^<-> ^^->) forwhich »a^.{e-(-^«r^A^->)}-min,{a.(z;a/->.A^->)}= 
Q^iz,a, ,A ) forwhich ^^^^ J^^^jg^^.^^,^^)}.^ 

This approach is implemented in the example embodiment by using a simulated 
annealing algorithm - a Monte Carlo approach for minimizing of multi-variable functions. 

The phase shift values thus determined are then utilised as phase shift 

values 5/ for the inverse scattering problem as part of the example embodiment (see equation 
(1)). The obtained dephased iV-channel design may be used as an input for the iterative 
optimisation procedure CTibodying the present invention. 

In the preferred embodiment, the strategy of difference minimisation is then extended 
further. The AT-central-channel-only minimisation set of dephasing angles ^, (I <i<N) is taken 
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and used to construct the initial sampling function S(z) = Q(z) e ^^'^ on a period 27i/M, Then the 
nontrivial Q(z) dependence is replaced by the asymptotic value N^^ and the resulting 
intermediate sampling function (with only phase sampling present) is decomposed in a Fourier 
series to find new Fourier coefficients at and a new set of the corresponding dephasing angles ^ 
s Arg(a/). This procedure leads to a change in partial higjh-amplitude grating coefficients. In 
iidditipn, n extra small higher order harmonic terms (corresponding to additional small 
anq>litude partial gratings) arise. Next we replace coefficients ai for N high-amplitude partial 
gratmgs with Ce , where a constant C is found firom the normaUsation condition 

t.M=NC' (4) 

and all coefficients a/ outside the band / e[l-/2/2,iV^ + n/2] being set to zero. These modified 
N+n amplitudes are used to construct a new "sampling*' function via inverse Fourier 
transformation. 

The procedure is repeated until a selected quality criterion has been reached. In the 
exanople embodiment, the iteration is stopped when a variation in the maximum "sampling" 
function amplitude between iterations falls below 10"^. 

The extended dephasing procedure for an infinite number of side channels n translates 
all nontrivial amplitude modulation of tiie sampling function into its phase ^z). As a price for 
doing that, about JON additional partial gratings (higher order harmonics with |a/|> 0.001) 
appear in the spectral characteristics, which means that phase dependence y^z) has a very fine 
structure. However, for practical applications the scale of the fine structure should always be 
larger than the size of the laser beam used for grating writing. An estimate for the upper limit 
for the maximum number N-^n of nonzero partial gratings can be given as 

125X1 

^ + ^= ' . (5) 

where Ao is the central wavelength [in microns], zq is the minimal laser beam size [in 
microns], AZ is the neighbouring channel spacing [in nanometers], and ho is the FBG average 
refiractive index. Estimate (5) was obtained using the assumption that about 8 grid points 
describe a period of the function cos(z) well enough. For the purpose of this description, the 
upper limit for the maximum number iSR-« given by the resolution limitation in the fabrication 
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of the grating [estimate (5)] will be referred to as un-limited in the context of the pieferred 
embodiment of the present invention, as it relates to a given experimental limitation. 

In the example embodiment, the nmnber of additional small amplitude partial gratings 
has been (finther) limited. It was found that selecting n^2N+2 channels gives good results in 
preferred embodiments of the invention. 

Accordingly, in the example embodiment the initial (pre-iteration) grating design 
fimction (compare Figure 1) is arrived at through a multi-chaimel grating design approach based 
on partial single grating dephasing, which may be further modified by an extended grating 
dephasing procedure allowing for a limited number of side channels in the spectral 
characteristics. The set of angles and amplitudes obtained by flie grating dq>hasing procedure is 
used for spectral dqihasing of partial channels to obtain a pre-iteration design. After this pre- 
iteration design is obtained, the iterative procedure embodying the present invention is q>plied. 

The implementation of the multi-channel grating design of the preferred embodiment in 
a gratmg structure requires grating writing apparatus with high spatial resolution to be utilised. 
Therefore, in a grating writing apparatus relying on photo induced refractive index changes, the 
apparatus preferably comprises a beam focusing means to reduce the size of the beam in the 
core of the photosensitive waveguide. 

Figure 7 shows an example experimental set up 50 for writing a multi-channel grating 
52 into an optical fibre 54. The experimental set up 50 comprises an interferometer 56 which 
includes a first acousto-optic modulation 58 being operated under an acousto-optic wave of a 
first fi^equency Qi, as indicated by arrow 14. An mcoming light beam 60 is incident on the first 
acousto-optic modulator 58 under a first order Bragg angle. The operating conditions of the 
acousto-optic modulator 58 are chos^ such that the modulator 58 is uifder driven, whereby 
approximately 50% of the incoming beam 60 is diffracted into a first order beam 62, and 50% 
passing through the acousto-optic modulator 58 as un-diffiacted beam 64. The un-diflBracted 
beam 64 is incident on a second acousto-optic modulator 66 of the interferometer 56 under a 
first order Bragg ang^ie, whereas the beam 62 is not. Accordingly, the beam 62 passing through 
the second acousto-optic modulator 66 without any significant loss. 

The second acousto-optic modulator 66 is operated under an acousto-optic wave of a 
fi'equency Q2, which propagates in a direction opposed the direction of the acousto-optic wave 
in the first modulator 58 as indicated by arrow 68. After the second acousto-optic modulator 66 
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the first order diffracted beam 70 and the beam 62 are frequency shifted in the same direction 
(e.g. higher frequency), but by different amounts i.e. Qi v Q2- 

The beams 62, 70 are then brougiht to interference utilising an optical lens 72, and the 
resulting interference pattern (at numeral 74) induces refractive index changes in the 
photosensitive optical fibre 54, whereby a refractive index profile, i.e. grating structure 52, is 
created in the optical fibre 54. 

In Figure 7, the optical fibre 54 is translated along the interferometer at a speed v, as 
indicated by arrow 74. 

It will be 25>preciated by a person skilled in the art that the experimental set up 50 shown 
in Figure 7 can be utilised to write a multi-channel grating structure of a multi-channel grating 
design embodying the present invention through suitable control of the first and second acousto- 
optic modulators 58, 66, in conjunction with a suitable control of the speed v at which the 
optical fibre 54 is translated along fiie interferometer 56 at any particular time. The higfh spatial 
resolution required to implement the multi-channel design of the preferred embodiment is 
achieved in the set up shown in Figure 7 by utilising optical lens 72, with the practical limit of 
the beam size in the focal plane preferably being of the order of the waveguide core size. 

It will be appreciated by a person skilled in the art that numerous variations and/or 
modifications may be made to the present invention as shown in the specific embodiments 
without departing firom the spirit of scope of the invention as broadly described. The present 
embodiments are, therefore, to be considered in all respects to be illustrative and not restrictive. 

For example, it will be appreciated that the initial grating design function utilised as the 
starting point of an improvement method embodying the present invention may be derived 
utilising different methods as the ones described for the example embodiment. Alternative 
methods can be found e.g. in Patent Co-operation Treaty (PCT) patent application no. 
PCT/AU02/00160 filed on 15 February 2002, entitled "Multi-channel grating design" assigned 
to the present applicant, PCT application no. PCT/AU02/00309 filed on 15 March 2002, 
entitled "Grating design" assigned to the present ^qpplicant, and/or Australian provisional patent 
application no. PS1131 filed on 15 March 2002, entitled ^Improved multi-channel grating 
design" in the name of the present applicant. 
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Furthennore, it will be appreciated tiiat the present invention is not limited to the type of 
multi-channel gratings of the example embodiments described. For example, the present 
invention can be equally applied to multi-channel gratings consisting of non-identical spectral 

r 

channels, and/or gratings consisting of groups of spectral chatmels. 

Furthermore, it will be appreciated that other numerical or ma&ematical approaches can 
be taken to inq)lement the iterative procedure embodying the present invention, e.g. comprises 
solving complex integral equations. 

Furthermore, multi-channel gratings can be created on the basis of the multi-channel 
grating design of the present invention using various known grating creation techniques, 
including one or more of the group of photo-induced refractive index variation in photo 
sensitive waveguide materials, etching techniques including etching techniques utilising a 
phasemask, and epitaxial techniques. 

In the claims that follow and in the smnmary of the invention, except where the context 
requires otherwise due to express language or necessary implication the word "comprising" is 
used in the sense of "including**, i.e. the features specified may be associated with further 
features in various embodiments of the invention. 
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Claims 



1. A metfiod of improving a grating design function describing a refractive index 



variation defining a multi-channel grating structure in a waveguide material, the improvement 



- modi^dng a first design function to generate a second design function having a reduced 
maximum amplitude compared with the first design function, 

- detCTnining a second response function associated with the second design function, 

- modifying the second response function to create a third response function having a 
desired functional spectral domain, and 

- determining a tiiird design function associated with the third response function, 

and iterating the method steps imtil the desired improvement is achieved, wherein the 
third design function of the previous iteration takes tiie place of the first design function of the 
next 

2. A method as claimed in claim 1, wherein the step of modifying the second 
response function comprises replacing the corresponding spectral domain of the second 
response function by the desired functional spectral domain. 

3. A method as claimed in claims 1 or 2, wherein the desired functional spectral 
domain comprises a corresponding spectral domain of a first response function associated with 
the first design function. 

4. A metiiod as claimed in any one of the preceding claims, wherein the method 
comprises the pre-iteration step of determining the first design function from a or the associated 
first response function. 

5. A method as claimed in any one of the preceding claims, wherein the 
determioing of the response functions from the design functions comprises solving a direct 
scattering problem, and the determining of the design functions from the response functions 
comprises solving an inverse scattering problem. 



being a reduced Trpayimiim refiractive index variation in the waveguide material along the 
grating stmcture while maintaining a desired functional spectral domain in a spectral response 
function associated with the design function, tiie method comprising the steps of: 
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6. A method as claimed in any one of the preceding claims, wherein the fliird 
response Amotion further has a desired response characteristic in at least one further spectral 
domain other than the fimctional spectral domain. 

7. A method as claimed in any one of the preceding claims, wherein the step of 
modifying the first design function comprises applying a normalisation process or an averaging 
process to the first design function. 

8. A method as claimed in claim 7, wherein the normalising process comprises 
replacing an amplitude function k{z) of the design function by a product of the square root of a 
constant A and a corresponding single-channel seeding amplitude fimction /r,(z), while 
maintaining a phase function of the design function. 

9. A method as claimed in claim 8, wherein the constant A is defined by the 

normalisation condition dz , where / is a length of the multi-channel grating 

0 0 

stmcture. 

10. A method as claimed in claim 7, wherein the averaging process comprises 
averaging over a sanntpling or quasi-sampling period of the design function. 

11. A method of fabricating a multi-channel grating stmcture, the method comprising 
the step of improving a grating design function describing a refiractive index variation defining 
the multi-channel grating stmcture in a waveguide material as claimed in any one of the 
preceding claims. 

12. A multi-channel grating stmcture fabricated utilising a method of fabrication as 
claimed in any one of flie preceding claims. 

Dated this 30th day of July 2002 



Redfern Optical Components Pty Ltd 
by its attorneys 
Freehills Carter Smith Beadle 
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